clc;
syms x;
y=x^2;
F=[x^2-2*x*y,y^2-2*x*y];
I=path_integral(F,[x,y],x,-1,1)